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Abstract 



We characterize a set of positive maps in matrix algebra of 4 x 4 complex matrices. Equivalently, 
we provide a subset of entanglement witnesses in C 4 ® C 4 parameterized by the rotation group SO(3) . 
Interestingly, these maps/witnesses define two intersecting convex cones in the 3-dimensional parameter 
space. The existence of two cones is related to the topological structure of the underlying orthogonal group. 
We perform detailed analysis of the corresponding geometric structure. 

1 Introduction 

One of the most important problems of quantum information theory is the characterization of mixed states of 
composed quantum systems P3H]- I n particular it is of primary importance to test whether a given quantum 
state exhibits quantum correlation, i.e. whether it is separable or entangled. 

The most general method of solving separability problem is the one based on the notion of positive maps or 
equivalently entanglement witnesses (EWs). A state p in Ha ® T~Lb is separable iff (1a <£> &)p > for all linear 
positive maps $ : B(Hb) B(Ha)- Recall, that a hermitian operator W € B(Ha £3> "Kb) is an entanglement 
witness [H [3] iff: i) it is not positively defined, i.e. W ^ 0, and ii) Tr(Wcr) > for all separable states a. 
Furthermore, a bipartite state p living in Ha <8> Hb is entangled iff there exists an EW W detecting this state, 
i.e. such that Tr (Wp) < 0. Due to the well known duality between maps and linear operators in B(Ha ®T~Lb) 
these two approaches are fully equivalent. Unfortunately, in spite of the considerable effort (see e.g. [5]-|18|L 
the structure of positive maps/entanglement witnesses is rather poorly understood. 

In this paper we analyze a class of positive maps $ : Mi(C) — > M^C) [M n (C) stands for an algebra n x n 
complex matrices] parameterized by the rotation group 50(3). This analysis extends our previous paper |14| 
where we discussed a class of maps $ : M^(C) — > Ms(C) parameterized by the commutative group SO (2). 
Our analysis shows that maps parameterized by SO (3) belong to two intersecting coaxial cones. We analyze 
the geometric structure of these convex. Interestingly, our construction recovers well known positive maps in 
Mi(C): reduction map and generalized Choi maps. We provide necessary and sufficient conditions for positivity 
and perform detailed analysis of (in)decomposability. Our discussion is illustrated by several geometric figures. 

It is hoped that our analysis sheds new light into the intricate structure of the convex cone of positive maps 
in matrix algebras. 

2 A class of positive maps in M n (C) 

Let us recall a construction of a class of positive maps in M n (C) introduced by Kossakowski in [10_ (for a slightly 
more general construction cf. [T^]). Let f a (a — 0, 1, . . . , n 2 — 1) denotes an orthonormal basis in M„(C), such 



that /o 



fk = fk, and 



Tr(/a//3) = S a p, a, (3 = 0, 1, . . . ,n 2 - 1 . 



(1) 



Note, that fk (k — 1, 2, . . . , n 2 — 1) are traceless, i.e. Trfk = . Now the positive unital map is defined as 
follows [TO] 

1 1 " 2_1 

$ fi (X) = -I n TrX + — - £ f a R a ^{ff}X) , (2) 

where i?^ is an arbitrary rotation matrix from 0(n 2 — 1). Note, that a dual map <I>^ defined by 

Tr[Y4> R (X)} = :Tv[^*(Y)-X] , 

reads 



I 1 " " 

(Y) = -I„Try + V f a R P MfpY) , (3) 

n n — 1 ^-^ 



a,/3=l 

and hence 

$ R (I n )=I n , $#(!„)= I„. (4) 
Note, that if i? corresponds to reflection in R™ _1 , i.e. i? Q) g = —S a p, one easily finds 

= — L^TrX - X) , (5) 

and hence it reproduced the reduction map in M„(C). 

Consider now a special class of maps corresponding to 

R = n® (-i» (n _i)) , (6) 

where 72. € 0(n — 1). In particular if 72. = — I„_i one reproduces reduction map. 

Let {ft} (k = 1, . . . , n 2 — 1) denote generalized Gell-Mann matrices defined as follows: let |1 ),..., \n) be an 
orthonormal basis in C" and define 



''' = ^- T (Ei fc )( fc i-^ + 1 )( / + 1 i) ■ ( ? ) 

vH' + L ) v fc= i 



for Z = l,...,n— 1, and 



uh = ^(|*)<i| + |0(*l), (8) 

vh = ^(\W\-\l){k\), (9) 

for k < I. It is clear that n 2 Hermitian matrices (/ , d/, itfc;, i^z) define a proper orthonormal basis in M n (C). 
One easily finds 

= E^'I^O'I' ( 10 ) 

*fl(|i>01) = -^yl*>01> Mi, (ii) 

where the matrix $jj reads as follows 

1 1 n_1 

= - + ^ E ^ (*l*|i> ■ (12) 
One shows [TT] that the matrix is doubly stochastic. 
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The corresponding entanglement witness is denned as follows 

W :=n(n- 1)(1 ® $r)P+ , 



(13) 



where P+ denotes maximally entangled state (we add the factor 'n(n — 1)' to simplify the form of W). One 
finds 



w= J2 \W\®Wij , 



where Wj, 



(j | for i j, and 



(14) 



W« := (n - l)$ R (\i)(i\) = (n - 1) £ 01 ■ 

3=1 

Example 1 If n = 3 one obtains the following formula for the 3x3 matrix <&ij 

^ I a b c 

*<j = 2 c ° 6 
y 6 c a 

where a, 6, c > are parameterized by the £0(2) rotation as follows 



2 , , 2 /, 1 \/3 . 

a = - ( 1 + cos a ) , o = — 1 cos a sin a 

3 v ; ' 3 I 2 2 



2/1 a/3 . 

c = - II — — cos a H — — sin a 



(15) 



(16) 



(17) 



TTiis class of maps was analyzed recently in [Hi. Note, that a + b + c = 2 and hence formulae define an 

ellipse on the bc-plane. Actually, this ellipse is defined by the following condition 



bc= {1-af 



(18) 



and for a < 1 it belongs to the boundary of a convex set of entanglement witnesses defined by the well known 
conditions HJjj 

i) 0<a<2, ii) a + b + c>2, Hi) a<l => be > (1 - a) 2 . (19) 

Moreover, it was shown f!5[ \161j that for a < 1 our family defines optimal entanglement witnesses. Interestingly, 
for a < 1 these witnesses are even exposed JJTjj. 



3 Entanglement witnesses in 



In this section we elaborate the construction of entanglement witnesses defined by ( 14 1 — ( 15 1 for n = 4. One 
finds for the corresponding $y matrix 



I ax bi ci dx ^ 

d 2 a 2 b 2 c 2 

c 3 d 3 a 3 b 3 

y 64 C4 0?4 04 J 



(20) 



where a,b,c,d > are parameterized by TZ € 50(3). Any such 7?. may be represented as follows 

TZ = -CpC~[S a — C a S~f C a CfjC^ — s a s^ C^Sp 



(21) 



where (a, /3,7) denote Euler angles and to simply notation we use s a := sin a and c Q := cos a. Unfortunately, 
the formulae for matrix elements <J>ij are quite involved (see the Appendix). Moreover, contrary to n = 3 
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the doubly stochastic matrix <frjj is no longer circulant. To simply our analysis we use the following simple 
observation: let Uki be a set of unitary matrices defined as follows 



(22) 



where uj = e 27 ™/ 4 and m © I denotes addition modulo 4. Let Pu = (I 4 ® Uki)P£ (I4 ® f/jy) and define 

W:=5^T>(Wflu)i^i . (23) 



It turns out [15] that if 14 7 is an entanglement witness then W is an entanglement witness as well. One finds 
that W is again defined by formulae (14 1— (15) with given by the following circulant matrix 



where 



/ a b c d \ 

d a b c 

c d a b 

\ b c d a J 



4 ^4 ^4 ^4 



(24) 



(25) 



read as follows 



3 + c a+7 (l + Cp) + Cp 



3 + -(s Q s 7 - CaC/j^ - 3c Q c 7 + 3cps 

ex. -2cp) 



2V3 



(3C-yS C 



2 2 

3c Q C / 3S 7 + CpC^Sa + C a Sy) + ^-^S^(2c 7 + C Q ) - -^SaS^ 



3 - -^(ZcaCpCy - 2s Q s 7 + cp) 



2 2 

^-^S^(C 7 - C a ) + -^ S P( S 1 + S a) 



1 

1 

4 



3 + -(s a S 7 — CotCpCry — iCaCj + 3cps a Sj — 2cp) -^{ScpC^Sc 

6 2y 3 



j 3Cqi tS^-y I Cf-y S Q£ I Cq C^Q S'-y ^ 



3V^ 



s,3(c 7 + 2c Q ) 



(26) 



(27) 



(28) 



(29) 



Interestingly, the EW corresponding to the reduction map does not belong to this class since reflection is not 
a proper rotation from SO(3). To include such case let us replace 1Z — > —1Z. It is clear that if 1Z is a proper 
rotation from £0(3) then —TZ € 0(3). Using the same arguments one obtains a new class of EWs W' with 
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(a, b, c, d) replaced by 



b' = 



d' 



3 - c Q+7 (l + cp) - eg 

3 ~^5o;S-y C(yCpc~y 3cQ.c-y I 3c^5Q/S-y " 2c^j ) " ^— (3c^ySQ I 
6 2v 3 



2 2 

+ 3c Q C,3S 7 + CpCjSa + C Q S 7 ) - — ^S^(2c 7 + C a ) + ^/g S a S /3 



1 2 2 

3 + ^faaCpOy - 2s Q s 7 + cp) + ^=s^(c 7 - c a ) - -^=s^(s 7 + s a ) + 



1 

1 



3 - ^(«qS 7 - c a c,gc 7 - 3c a c 7 + 3cps a s^ - 2cg) + ^=(3cpc 7 s a + 

2 2 

+ 3c Q s 7 + c 7 s Q + c a cps y ) + — ^=s ( 3(c 7 + 2c Q ) + — ^=s 7 s^ 



Formulae for (a, 6, c, d) and (a', 6', c', d') provide an analog of much simpler relations ( 17) for n = 3 



(30) 



(31) 



(32) 



(33) 



4 Geometry of cones 



Now comes a natural question: what is the geometric representation of the above relations? For n = 3 formulae 



(17) give rise to an ellipse in the 6c-plane (cf. [3]). Interestingly for n — 4 an elegant geometric picture arises as 
well. Actually, it was the original motivation for this paper. Numerical analysis shows the following picture in 



(6, c, d) coordinates (recall, that a + b + c + d = 3): formulae ( 26 1— ( 29 ) and ( 30 1 ( 33 1 give rise to two intersecting 
cones (see Fig. (a)) The above cones are described by the following equations: Cone I 



(1/2,1,1/2)+(1,-1,1)t 





(a) 

Figure 1: Cone I [blue] is bounded by a plane b = —d + 2 and cone II [red] by a plane b = —d + 1. Both 
cones are coaxial with axis described by the line (|, 1, \) + (1, —1, l)t. Moreover, the special points are marked: 
associated with generalized Choi maps (i), (ii), (iv) and connected with reduction map (iii). 

(b - if + (2c - 3) 2 + (d- 2) 2 + 46c + Acd - 2bd = 9 , (34) 

and Cone II 

(b - l) 2 + (2c - 3) 2 + (d- l) 2 + 46c + Acd - 2bd = 6 , (35) 
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with the constraints that — 6 + 1 < e? < — 6 + 2. The vertices of these cones are located at b = \, c = 1, d = h 
and 6 = 1, c = \, d = 1 for cones I and II, respectively. They intersect along an ellipse in the plane d = — 6+ ¥. 



Let us analyze the intersection of the Cone I defined by (34) with the plane 6 + d = 2. One finds 



(d-lf + {c- l -f = \ . (36) 



Similarly the intersection of the Cone II defined by ( 35 ) with the plane 6 + d = 1 gives 



(c-l) 2 + (d-^) 2 = \ • (37) 



Taking into account a+b+c+d=3 one finds 



a = 1 — c, 

6 = 1± Vc(l-c), (38) 
d = lT#-c), 



with c e [0, 1] for the ellipse I on 6 + d = 2 plane and 



a = 1 ± Vd(l-d), 

6 = l-d, (39) 
c 



= 1 + - d) 

with d € [0, 1] for the ellipse II on 6 + d = 1 plane. Interestingly, ellipse I [blue] satisfies 

a + c = 1 , b + d = 2 , (40) 



whereas the ellipse II [red] satisfies 



Formulae ( 39 1 imply 



Similarly, formulae ( 38 1 imply 



a + c = 2, b + d=l. (41) 
W = (1 - a) 2 . (42) 
ac = (1 - 6) 2 . (43) 



These two equations provide an analog of the well known condition ( 18 I for n — 3 



Let us observe that these two ellipses contain already known positive maps: 

• $[0, 1, 1, 1] - reduction map [point (iii) in Fig. 1], 

• $[1, 1, 1, 0] - generalized Choi map [point (i) in Fig. 1], 

• $[1,0, 1, 1] - generalized Choi map [point (ii) in Fig. 1], 

• $[1,1,0,1] - [point (iv) in Fig. 1[. 

Note, however, that another generalized Choi map $[2, 1, 0, 0] does not belong to our class. 

5 (In)decomposability 

In this section we analyze the issue of indecomposability of $[a, 6, c, d]. 
Theorem 1 $[a,6, c, d] is decomposable if and only if b = d. 
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Proof: let us consider a state given by (unnormalized) density matrix 



p s =J2 []**><**! + + \i,i + 2)(i,i + 2\ + s- 1 \i,i + 3)(i,i + 3|] + , (44) 

i=l i^j 

with e > 0. One easily checks that p e is PPT. One finds 

Tr(W[a, 6, c, d]p e ) = -12 + 4(a + c) + Ae^d + Aeb = 4[de~ 1 + be - (b + d)] = ^[be 2 - (b + d)e + d] , (45) 

and hence Ti(W[a,b,c,d]p E ) < iff there exists e > such that be 2 — (b + d)e + d < 0. The corresponding 
discriminant reads 

A= {b + d) 2 -4bd= (b~d) 2 , 
and hence be 2 — (6 + d)e + d < if e £ (e_,e + ) with 

b + d±\b-d\ 



2b 



Note, that e + > e_ if and only if b ^ d. This way we have proved that if b 7^ d then $[a, 6, c, d] is indecomposable. 
Now we show that if b = d, then $[a, 6, c, d] is decomposable. We find the corresponding decomposition 



W[a, 6, c, d] — P[a, 6, c, d] + Q[a, b, c, d] T , 

where, P[a,b,c,d] and Q[a, b, c, cf| are positive operators. One has 
4 

P[a, b, c, d] = [o|»> (»| - (1 - 6) [ |*») <* + 1, * + 1| + I it) (t + 3, i + 3| ] - (1 - c)|m) (i + 2, i + 2| 



(46) 



(47) 



and 



Q[a,6,c,d] = ^^|i,i + l)(i,i + l| + c|i,i + 2)(i,i + 2|+6|i,i + 3)(i,i + 3| 



b(|i,i + l)<i + l,i| + \i,i + 3){i + 3,i\)-c\i,i + 2){i + 2,i\ 



(48) 



A = 



(49) 



It is clear that Q[a, b, c, d] > 0. Now, to prove that P[a, b, c, d] > one needs to show that the following circulant 
matrix 

/ a 6-1 c-1 6-1 \ 
6-1 a 6-1 c-1 
c-1 6-1 a 6-1 
\ 6-1 c- 1 6-1 a J 

is positive. One finds for the eigenvalues of A: {0, 4(1 — 6), 2(2 — 6 — c), 2(2 — 6 — c)}. Note, that 6 + d < 2 and 
hence 6 = d < 1. Moreover, 2 — 6 — c=(a + c) — 1 > 0. Hence A > 0. 

□ 



6 Structural Physical Approximation 

For any entanglement witness W in T-La ®'Hb such that Tr l-F = 1 one defines its structural physical approxi- 
mation (SPA) 

W{ P ) = (1- P )W + -^—l A ®l Bl (50) 

with p > p* , where p* is the smallest value of p such that W(p) > 0. Thus, SPA defines a legitimate quantum 
state in T-La^T^b- The conjecture of Korbicz et al. [5D] (see also [2T]) states that, if W is an optimal EW, then 
its SPA defines separable state. It was supported by several examples (see e.g. |22|). 



In a recent paper [14] it was conjectured that all entanglement witnesses W[a, 6, c] with a, 6, c satisfying ( 18 ) 
and a < 1 are optimal. Actually, this conjecture was proved by Ha and Kye [T5] (see also [IS]). It was shown 
[Tl"] that W[a, 6, c] support SPA conjecture |5D]. Now we prove the following 
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Figure 2: The intersection of the plane b = d and of the blue and red cones represents decomposable EWs. 



Proposition 1 If a, 6, c, d satisfy or (38) then the structural physical approximation to W[a, b, c, d] provides 
a separable state in C 4 <8> C 4 . 



Proof: Let us consider SPA for our class 

w(p) = (i- P )w + ^i 4 ®]u . 

Now, W(p) > for p> p* , where the critical value p* is given by 

4(a - 3) 



P 



It turns out that W(jf) can be represented as 

1 



W(j3* 



4[3 + 4(o-3)] 



3 + 4(a-3)' 

(CTi2 + 0"i3 + (Ti4 + (723 + 024 + 034 + O d) 



where 
and 



(51) 
(52) 

(53) 
(54) 



4 

^ = 51 (C2b + c + d-l)\i,i + l)(i,i + l\ + (2c + b + d-l)\i,i + 2)(i,i + 2\ + (2d + b + c-l)\i,i + 3)(i,i + 3\y (55) 

i=l 

Due to the fact that try are PPT and supported on C 2 ® C 2 , they are separable. Moreover, ad is separable, 
whenever it defines a legitimate quantum state, that is, when 



2b + c + d > 1, 
2c + b + d > 1, 
2d + b + c > 1. 



(56) 



It is straightforward to show that both conditions ( 39 1 and ( 38 ) imply ( 56 1 which ends the proof. 
Interestingly, the above three 2-dimensional planes: 

2b + c + d = 1 , 2c + b + d = 1 , 2<i + b + c = 1 , 

intersect at b = c = d = 1/4. 



□ 
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7 Conclusions 



We analyzed a class of positive maps in M^C) (or equivalently class of EWs in £>(C 4 <g> C 4 )). This construction 
generalizes analysis provided in [14 for n = 3 by replacing commutative group SO (2) by the noncommutative 
rotation group 50(3). We formulated necessary and sufficient conditions for positivity of $[a,b,c,d] and de- 
scribed the geometric structure of the convex set formed by these maps. Interestingly, there are two proper cones 
in the 3-dimensional space parameterized by (b, c, d) . It was shown that for b ^ d all maps are indecomposable 



and hence can be used to detect PPT entangled states. Moreover, maps satisfying (56) support SPA conjecture 



[20] . We provided two natural 1-parameter subclasses of maps (cf. formulae (39) and (38 1 ) - two ellipses in the 
parameter space - which are direct generalizations of 1-parameter class of maps analyzed in |14| . 

In a forthcoming paper we are going to analyze further properties of maps $ [a, 6, c, d] . In particular optimality 
and atomicity. Finally, it would be interesting to provide the analogous construction for arbitrary n. 

Appendix 

a 1 = l + \ Rll + ^=R 12 + ^=R 13 + ^=R 21 + \r 22 + ^=R 23 + ^R 31 + ^R 32 + ^i? 33 , 

a 3 = 3 - - -^R 32 + ±R 33 + 2 -R 22 -^R 23 
3 4 3V2 3 12 33 + 3 K22 3^2 23 ' 

3 9 

° 4 ~ 4 + 12 R33 ' 

61 = 1 " 2 fll1 + sTi^ 12 + 27i i?13 " 271 i?21 + e^ 22 + ^ 23 - + ^ R32 + i^ 33 ' 

b 2 = 3 - + ^R 12 - -^R 13 - l -R 22 + -^R 23 ~ -^R 32 + ±R 33 

2 4 + ^3 2V6 3 22 6V2 23 3^2 32 12^ 33 ' 

63 = 4-^33 + ^23, 

64 = i " 2 7i i?31 " 27i i?32 _ T^ 33 ' 

* = - - — fl U + — ilia ~ -R22 + — i? 23 - — R 3 2 + -i? 33 
1 4 V3 12 2^6 " 3 22 6V2 23 3^2 32 12^ 33 ' 

, - 3 1 ff 4- 3 R 3 7? 

C2 " i " 271 23 + 27! 13 " ^ 33 ' 

3 4 + 2V6 31 6^2 32 12 33 V3 21 3^ 22 3^2 23 ' 

, - 3 + 3 R 1 R 3 R 
C4 " 4 + 271 31 " 271 32 " 12^ 33 ' 

dl = 4 ~ 27^" " 2^" " 4 fl33 ' 

d2 = 4 + 271 i?21 + e^ 22 + eT^ 23 ~ \ Rl1 ~ ^ Rl2 ~ 2VE Rl3 + 2VE Ral + eT^ 32 + u R ™ ' 
ds = 3 - — R*i + — + -R 33 + -R 2 i - -R 22 - —R 23 

3 4 2V6 6V2 32 12 33 + VS 3 22 3V2 23 ' 

di = 4 + T^ 32 " 12^ 33 • 
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